The Cumulative Bakground of High-Energy Neutrinos from Starburst Galaxies 
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We show that starburst galaxies convert efficiently cosmic-rays into pions, which in turn decay 
into high-energy neutrinos and photons. The cumulative background of GeV neutrinos is E 2 ^> v « 
hT 7 GeV cm 2 s 1 sr \ and its extrapolation to higher neutrino energies depends on the energy 
spectrum of the injected cosmic-rays, with El<& v cx EZ °- ls±0 - 1 up to E v ~ 0.3 PeV and possibly 
higher neutrino energies. This flux, which constitutes a lower limit to the high energy extra-Galactic 
neutrino flux, is potentially detectable by forthcoming km-scale neutrino telescopes. 

PACS numbers: 95.85.Ry, 98.54.Ep, 98.70.Sa, 14.60.Pq 



Large area, high-energy neutrino telescopes are be- 
ingconstructed to detect extra-Galactic neutrino sources 
[l|, |2j • It is widely believed that high-energy neutrinos are 
produced at cosmological distances based on the fact that 
the cosmic-ray energy spectrum extends to > 10 20 eV and 
is most likely dominated above ~3x 10 18 eV by extra- 
Galactic protons 0- High-energy neutrinos are then 
likely to be emitted by the same sources that produce 
the high-energy protons, through the decay of charged 
pions generated by interactions of these protons with 
the radiation field (or nucleons) within these sources . 
Gamma-ray bursts (GRBs,@) and jets of active galactic 
nuclei (AGN, 0) have been suggested as possible sources 
of high-energy neutrinos that are associated with high- 
energy cosmic-rays. 

If the source size is not much larger than the mean- 
free-path of ultra high-energy cosmic-rays (UHECRs) for 
pion production, as is the case for both AGN jets and 
GRBs, then UHECR observations set an upper bound of 
£ 2 $„ < El<S>y B = 5 x 10~ 8 GeV cm" 2 s" 1 sr" 1 on the 
flux of high-energy neutrinos produced by pion decay [6( . 
This so-called Waxman-Bahcall (WB) upper bound, im- 
plies that km (Gton) scale neutrino detectors are required 
for detecting the high energy neutrino flux produced by 
sources of UHECRs in the energy range of ~ 1 TeV to 
~ 1 PeV .and larger detectors are required at yet higher 
energies p. 

The extra-Galactic neutrino number flux per unit en- 
ergy is expected to be close to $^ VB in the energy range 
of 10 19 -T0 20 eV .7], since protons originating at cosmo- 
logical distances with > 5 x 10 19 eV lose all their energy 
to pion production through their interaction with cosmic 
microwave background (CMB) photons [Hj. However, 
the existence of a neutrino flux at the level of for 
lower energies is not guaranteed. Protons in this energy 
range may lose only a small fraction of their energy to 
pion production, making a neutrino flux much lower than 

$WB 

In this Letter we show that radio observations of star- 
burst galaxies imply a lower limit on their cumulative 



extra-Galactic neutrino background flux <&^ B which is of 
order at the energy range of ~ 1 GeV to ~ 300 TeV 
(0.3 PeV), and possibly extending also to higher ener- 
gies. Previous discussions |9( have used a different model- 
dependent approach to estimate the neutrino luminosi- 
ties of individual galaxies and did not integrate the back- 
ground neutrino flux over cosmic history. The main un- 
certainty in our estimate involves the spectral index of 
the injected cosmic-rays. As we show next, starburst 
galaxies are " hidden" cosmic-ray sources in the sense that 
they dissipate the cosmic-ray energy produced within 
them, and so their neutrino background could in prin- 
ciple exceed the WB bound. Nevertheless, the neutrino 
flux we derive is comparable to the WB bound, because 
the energy production rate of cosmic rays in starbursts 
is close to the cosmic production rate of UHECRs. As 
discussed below, this is not necessarily a coincidence. 

The detection of multi-TeV neutrinos by upcoming 
km-scale detectors would provide unique constraints on 
models for cosmic-ray generation, as well as a unique han- 
dle on the physical properties of the interstellar medium 
in starburst galaxies. The expected detection rate may 
allow km scale detectors to study the flavor ratio of as- 
trophysical neutrinos, thus probing neutrino oscillations 
and physics across unprecedented scales of energies and 
lengths 

Starbursts as neutrino factories. A substantial frac- 
tion of the cosmic star formation activity at redshift z ~ 2 
[I2I Il3j | occurs in transient starburst episodes. These 
episodes are often triggered by galaxy merger events, 
which channel fresh gas towards the center of the merger 
remnant 0, ^| . The characteristic scale over which the 
gas is observed to be concentrated, (£ ~ hundreds of pc) 
and the typical gas velocities, (v ~few hundred km s _1 ) 
imply a dynamical time of a few million years, compa- 
rable to the lifetime of massive stars. Core-collapse su- 
pernovae (SNe) are therefore expected to enrich the gas 
with relativistic protons and electrons, i.e. cosmic-rays, 
which are accelerated in the collisionless shocks produced 
by these explosions. Additional relativistic particles may 
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be injected by an outflow from a central supermassive 
black hole, in cases where quasar activity accompanies 
the starburst phase 0, . Synchrotron radio emission 
is routinely observed in starburst galaxies, confirming the 
presence of relativistic electrons within them 01 • 

Relativistic protons, injected along with the electrons 
into the starburst interstellar medium, would lose energy 
primarily through pion production by inelastic collisions 
with interstellar nucleons. The decay of charged pions, 
7T+ — * fi + + — > e + + v e + + and 7r~ — > fj,~ + 9^ — > 
e~ + v e + + Up, would then convert part of the proton 
energy to neutrinos. In what follows we estimate the 
neutrino flux by showing that: (i) protons lose essentially 
all their energy to pion production, and (ii) a lower limit 
to the energy loss rate of protons may be directly derived 
from the synchrotron radio flux of the secondary e^. 

Protons would lose all their energy to pion produc- 
tion provided that the energy loss time is shorter than 
both the starburst lifetime and the magnetic confine- 
ment time within the starburst gas. In the energy range 
of interest, the inelastic nuclear collision cross section is 
w 50 mb, with inelasticity of w 0.5 The energy 

loss time, ti oss (0.5ncr pp c) _:L where n is the interstel- 
lar nucleon density, would be shorter than the starburst 
lifetime which is at least a dynamical time, ~ (2£/v), as 
long as 

Sgas > Scrit = ^=3X 10- 2 /3_ 3 g Cm- 2 . (1) 
CTpp 

Here £ gas ~ rn p n£ is the surface mass density of the gas 
and (3 = v/c = /3_3(u/300 km s -1 ). As it turns out, 
the critical surface density, S cr ;t, is comparable to the 
minimum £ gas in known starburst galaxies |l9j |. 

The ratio of confinement time, r con f, and loss time, 
ti oss , is less straightforward to estimate, since magnetic 
confinement of cosmic-rays is not well understood. For 
starburst galaxies £ gas pa 10 9 M Q /kpc 2 = 0.2g cm -2 , 
which for a disk height of ~ 100 pc |20| implies 
Tioss ~ 10 5 yr. The confinement time of 10 GeV pro- 
tons in our Galaxy is r con f ~ 10 7 yr |2l|. The to- 
tal gas column density traversed by protons of energy 
< 1 TeV before they escape the Galaxy is T^ on { « 
9(£/10GeV)- s g cm- 2 with 0.5 < s < 0.6 22], sug- 
gesting r conf ~ 10 7 (-B/10GeV)-°- 6 yr for the most 
commonly-used value of s = 0.6. If the confinement time 
in starburst galaxies was similar to that of the Galaxy, 
then the comparison of ti oss with r con f would have im- 
plied that protons with E < 30 TeV lose all their en- 
ergy prior to their escape from starbursts. However, the 
magnetic field in starbursts is much larger than in the 
Galaxy. Thompson et al. 0] have found that the mag- 
netic field strength within starburst galaxies scales lin- 
early with Sgas, and has an amplitude that is ~ 100 
times larger than in the Galaxy for £ gas = 0.2g cm- 2 . 
Assuming that the confinement time depends on energy 
only through the proton's Larmor radius (oc E/B), this 



suggests that protons with E < 3 x 10 3 TeV lose all their 
energy prior to their escape from typical starburst galax- 
ies. Moreover, the neutrino flux is expected to be domi- 
nated by starbursts at z > 1, for which the typical surface 
density should be even higher than in local starbursts. It 
is therefore reasonable to assume that most of the en- 
ergy injected into starburst galaxies in E < 3 x 10 3 TeV 
protons is converted to pions. 

We now proceed to calibrate the luminosities of star- 
burst galaxies in high-energy neutrinos based on their 
observed synchrotron luminosities. As recently shown 
by Thompson et al. the synchrotron cooling time 

of relativistic electrons radiating at a frequency of few 
GHz is typically much shorter than the starburst lifetime. 
The synchrotron luminosity reflects, therefore, the energy 
production rate in these electrons. The production rate 
of electrons with energy E e , E^dN e /dE e , is related to 
the synchrotron radio luminosity L v per unit frequency v 
produced by these electrons, as E^dN / dE e « 2vL v . The 
factor of 2 originates from the fact that the synchrotron 
frequency is proportional to the square of the electron 
Lorentz factor, v oc -E 2 , so that the energy of electrons in 
a given decade of E e is spread over two decades of pho- 
ton frequency v. Denoting the ratio of injected power 
in protons and electron at a fixed particle energy as 
r/p/e, the luminosity in and per logarithmic energy 
bin E v is given by E u dL/dE u = (l/3)ri p/e E^dN /dE e = 
(2/3)rj p / e b'L l/ . The factor of 1/3 comes from the fact that 
w 2/3 of the proton energy is carried by charged pions 
(~ 1/3 by neutral pions), and ps 1/2 of the pion energy 
is carried, after it decays, by muon neutrinos. Cosmic- 
ray data in the Milky- Way galaxy indicates r\ p j e ~ 50 
(e.g. H3)- However, we expect a lower ratio in starburst 
galaxies. The decay of the charged pions produces sec- 
ondary e ± which carry w 2/3 x 1/4= 1/6 of the proton 
energy (since the carry ps 1/4 of the decaying tt^ en- 
ergy). This implies rj p / e rj 6 and a neutrino luminosity 
E v dLj dE v = AvL u . We conservatively ignore the ion- 
ization loss by the electrons whose inclusion would 
only increase our predicted neutrino flux. 

In the Milky way, electrons radiating at few GHz have 
energies E e ~ 10 GeV. In starburst galaxies the magnetic 
field is typically ~ 100 times larger, implying that the ra- 
dio luminosity observed is produced by E e ~ 1 GeV elec- 
trons. Since the electrons and neutrinos produced in pion 
decays carry similar energy per particle, the radio lumi- 
nosity provides a direct estimate of the GeV neutrino lu- 
minosity, EydL/dE v {E y — lGeV) = 4zaL„|„^ighz- The 
local 1.4 GHz energy production rate per unit volume is 
v{dL v /dV)\ y=1AG Kz ps 10 43 erg Mp C - 3 yr" 1 17j. Most 
of the stars in the Universe formed at redshifts z ~ 2-4 
[2^ | inside starburst galaxies 0, 0] . The fact that the 
starburst mode of star formation dominated at z > 2 
is also supported by the prominence of quasar activity 
around the same redshift, as both processes are believed 
to result from mergers of gas-rich galaxies. We extrap- 
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olate the local 1.4 GHz energy production rate per unit 
volume (of which a dominant fraction is produced in qui- 
escent spiral galaxies) to the redshifts where most of the 
stars had formed through the starburst mode, based on 
the observed redshift evolution of the cosmic star forma- 
tion rate [24|, and calculate the resulting neutrino back- 
ground. The cumulative GeV neutrino background from 
starburst galaxies is then 



E 2 ME V = lGcV) « — (t H [4K^/W)Ui.4GH Z 



= l(T Y Co.5 GeV cm -2 s" 1 sr" 1 



(2) 



Here, tn is the age of the Universe, and the factor 
( = 10°' 5 (o.5 incorporates a correction due to redshift 
evolution of the star formation rate relative to its present- 
day value. The value of Co. 5 ~ 1 applies to activity that 
traces the cosmic star formation history [jj. Note that 
flavor oscillations would convert the pion decay flavor ra- 
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FIG. 1: The shaded region brackets the range of plausible 
choices for the spectrum of the neutrino background. Its up- 
per boundary is obtained for a power-law index p — 2 of 
the injected cosmic-rays, and its lower boundary corresponds 
to p = 2.25 for E v < 10 14 ' 5 eV. The solid green line corre- 
sponds to the likely value p = 2.15 (see text). Other lines: the 
WB upper bound on the high energy muon neutrino intensity 
from optically-thin sources; the neutrino intensity expected 
from interaction with CMB photons (GZK); the atmospheric 
neutrino background; experimental upper bounds of optical 
Cerenkov experiments (BAIKAL and AMANDA |3^): 
and the expected sensitivity of 0.1 km 2 and 1 km 2 optical 
Cerenkov detectors 0. 

Equation J3J) provides an estimate of the GeV neu- 
trino background. The extrapolation of this background 
to higher neutrino energies depends on the energy spec- 
trum of the high energy protons. If the proton energy dis- 
tribution follows a power-law, dN/dE cx E~ p , then the 



neutrino spectrum would be, E^$> v oc E 2 ~ p . The energy 
distribution of cosmic-ray protons measured on Earth fol- 
lows a power-law dN/dE oc E~ 2J5 up to the "knee" in 
the cosmic-ray spectrum at a few times 10 15 eV [13, S3- 
(The proton spectrum becomes steeper, i.e. softer, at 
higher energies 0.) Given the energy dependence of the 
confinement time, oc E~ s |22| . this implies a produc- 
tion spectrum dN/dE oc E~ p with p = 2.75 - s » 2.15. 
This power-law index is close to, but somewhat higher 
than, the theoretical value p — 2, which implies equal 
energy per logarithmic particle energy bin, obtained for 
Fermi acceleration in strong shocks under the test par- 
ticle approximation [2^|. We note that the cosmic-ray 
spectrum observed on Earth may not be representative 
of the cosmic-ray distribution in the Galaxy in general. 
The inferred excess relative to model predictions of the 
> 1 GeV photon flux from the inner Galaxy, implies that 
the cosmic-rays are generated with a spectral index p 
smaller than the value p — 2.15 inferred from the local 
cosmic-ray distribution, and possibly that the spectral 
index of cosmic-rays in the inner Galaxy is smaller than 
the local one |27J . The spectrum of electrons accelerated 
in SNe is inferred to be a power law with spectral index 
p = 2.1 ± 0.1 over a wide range energies, ~ 1 GeV to 
~ 10 TeV, based on radio, X-ray and TeV observations 
(e.g. 0). 

For a steeply falling proton spectrum such as dN/dE ~ 
E~ 2 , the production of neutrinos of energy E v is domi- 
nated by protons of energy E « 2QE V |l8|, so that the 
cosmic-ray "knee" corresponds to E v ~ 0.1 PeV. In anal- 
ogy with the Galactic injection parameters of cosmic- 
rays, we expect the neutrino background to scale as 



£ 2 $f « 10- 7 (^/lGeV)-° 15±01 GeV cm- 2 s" 1 sr^S) 



up to ~ 0.1 PeV. In fact, the "knee" in the proton spec- 
trum for starburst galaxies may occur at an energy higher 
than in the Galaxy. The steepening (softening) of the 
proton spectrum at the knee may be either due to a 
steeper proton production spectrum at higher energies, or 
a faster decline with energy for the proton confinement 
time. Since both the acceleration of protons and their 
confinement depend on the magnetic field, we expect the 
"knee" to shift to a higher energy in starbursts, where the 
magnetic field is much stronger than the Galactic value. 
The predicted neutrino intensity is shown as a solid line 
in Fig. ^ The shaded region illustrating the range of 
uncertainty in the predicted neutrino background. This 
range is bounded from above by the intensity obtained 
for p = 2, corresponding to equal proton energy per log- 
arithmic bin, and from below by the intensity obtained 
for p — 2.25, corresponding to the lower value of the 
confinement time spectral index, s = 0.5. 

The extension of the neutrino spectrum to energies 
E u > 1 PeV is highly uncertain. If the steepening of the 
proton spectrum at the knee is due to a rapid decrease 
in the proton confinement time within the Galaxy rather 
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than a change in the production spectrum, then the neu- 
trino background may follow Eq. Q above 1 PeV, pro- 
vided that the magnetic confinement time in starbursts 
does not decline rapidly at these energies. This exten- 
sion is shown as the dashed line in Fig. ^ If, on the other 
hand, the steepening at the knee is due to a steepening in 
the production spectrum, then the neutrino background 
may decrease rapidly following the lower boundary of the 
shaded region. Note that although transient sources such 
as GRBs or AGN may not contribute at the present-time 
to the observed Galactic cosmic-rays, they could domi- 
nate the cosmic integral of the neutrino emission by star- 
burst galaxies 31]. Such a contribution would be par- 
ticularly important at E v > 0.1 PeV if the same sources 
produce the UHECRs. 

Gamma-ray Production. The neutrino emission from 
7T^~ decay will be accompanied by a comparable flux of 
7-rays from ir° decay with energies E 1 > lGeV. Since 
the protons carry more energy than the electrons and 
the electrons lose only a minor fraction of their energy to 
7-rays (through inverse-Compton scattering of soft pho- 
tons |l9|), the 7-rays provide a robust measure of the 
pion production rate by protons. This allows a sanity 
check on our predicted neutrino fluxes. Our expected 
7-ray flux from the nearest starburst galaxies must be 
smaller than existing upper-limits. For example, analy- 
sis of HESS and EGRET data for NGC253 provided an 
upper limit of ~ 1.9 x 10 -12 erg cm -2 s _1 on its 7-ray 
flux as a point source at E 1 > 1 TeV 0, ■ Based on 
its 1.4GHz flux of 10 -13 erg cm -2 s _1 , we predict a 7- 
ray flux of < 10~ 12 erg cm' 2 s -1 , well below the existing 
upper limit. We have repeated this calculation for other 
starburst galaxies with known upper limits || , and found 
similar results. Our expected 7-ray fluxes are consistent 
with previous models for individual sources which 
were based on other considerations. The forthcoming 
GLAST [34| mission will have the sensitivity to detect 
the predicted 7-ray fluxes from nearby starbursts and 
revalidate our expectation for the neutrino background. 

Discussion. The neutrino intensity predicted by 
Eq. © at 100 TeV is w 2 x 10" 8±0 5 Gev cm" 2 s" 1 sr" 1 , 
implying a detection rate at E v > 100 TeV of w io 15±0 - 5 
events per yr in a 1 km 2 detector. The neutrino back- 
ground is therefore potentially detectable by forthcoming 
km-scale telescopes. If detected, the signal would serve 
as a highly-effective probe of cosmic-ray production in 
star forming environments. 

As already mentioned, our inferred background of 
high-energy neutrinos is similar to $ WB since the 
present-day energy production rate of cosmic-ray protons 
by starbursts, E 2 dn/dE = 12v{dL v /dV)\u=x.^GOz ~ 
10 44 erg Mpc~ 3 yr" 1 , is similar to the energy produc- 
tion rate of UHECRs |33j. While the neutrino back- 
ground originates primarily in starburst galaxies with 
£ gas ^ ^crit, extra-Galactic cosmic-rays would originate 
in galaxies with S gas < E cr it, where they do not en- 



counter substantial loses on their way out. Since the two 
classes of galaxies are responsible for making comparable 
fractions (up to a factor of ~ 2) of the stellar mass reser- 
voir in the local universe [lii ]. the coincidence between 
and $wb 

may simply reflect the matching condi- 
tion at E gas ~ S cr it that known galactic environments 
account for the extragalactic flux of UHECRs. Thus, as 
long as the sources of UHECRs (such as GRBs or AGN) 
exist in starburst galaxies, they would produce a mea- 
surable background of neutrinos with energies E v > 10 
eV (near the upper envelope of the shaded region in Fig. 
1), that extends well beyond the atomspheric neutrino 
background. 

Our results do not depend on whether the cosmic- 
rays originate from supernovae or the supermassive black 
holes in galactic nuclei. While supernova-produced 7- 
rays will be spread throughout the star-forming volume 
of these nuclei, a black hole origin would bias the 7-rays 
to be centrally concentrated. GLAST will be able to dis- 
tinguish between these possibilities by imaging nearby 
starburst galaxies in 7-rays. 
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